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Summary 


Current spray materials utilized by the USDA Forest Service are non-Newtonian 
in nature; a better understanding of their atomization behavior will improve their 
efficiency, efficacy, and total accountantacy and environmental fate. Thus, this study 
first reviews the state of Newtonian and non-Newtonian spray material characteristics, 
and attempts data compression with available wind tunnel data. Starting with first 
principles of dimensional analysis, we develop consistent nondimensional parameters that 
group into well-accepted numbers, such as the Reynolds and Weber numbers. We 
correlate the water atomization test runs from the USDA Forest Service drop size 
distribution data base with simple expressions that recover the details of the spray. Then 
we speculate upon the form of the non-Newtonian terms for the atomization of 
viscoelastic fluids. Finally, we suggest that a systematic laboratory benchtop 
examination of applicable non-Newtonian materials be undertaken, followed by a 
consistent set of wind tunnel tests to encompass anticipated spray materials and field 
conditions. 


1. Introduction 


_ The drop size distribution of spray material atomized by nozzles influences the 
magnitude of evaporation, spray deposition, drift, and application effectiveness. Several 
factors affecting the atomization process and the resulting drop size distribution include 
nozzle type, hydraulic line pressure, flow rate, shear across the nozzle orifice, and 
specific properties of the formulation. The details of the resulting drop size distribution 
are critical to forest and agricultural applications, where specific levels of spray material 
must be deposited to achieve success and avoid environmental insult. 


In an effort to build a data base of typical formulations and aerial application 
conditions, the USDA Forest Service, and other agencies and companies, have contracted 
over the last several years for wind tunnel tests to determine drop size distributions of 
pesticides and simulant spray material atomized by nozzles. The purpose of these studies 
is to determine the atomization of tank mixes as influenced by hydraulic line pressure, 
flow rate, air velocity and shear across the atomizer, components including adjuvants 
(chemical, physical and biological) in the tank mix, viscosity, specific gravity, surface 
tension, and other atmospheric conditions. The intent is to understand the behavior of 
these spray materials not only to support registration and re-registration with the U. S. 
Environmental Protection Agency, but also to reduce the risk of material drift. 


Current spray materials utilized by the USDA Forest Service are non-Newtonian 
in nature; a better understanding of their atomization behavior will improve their 
efficiency, efficacy, and total accountantacy and environmental fate. Thus, the driftable 
part of the atomization must be known; however, this step at present would require a 
wind tunnel test. Therefore, it becomes prudent to seek approaches that might reduce the 
need for such wind tunnel testing. One such approach is to develop a limited data base of 
wind tunnel tests, correlate these data, and then use them to extend predictions to 
conditions not tested. 


A data base has already been compiled (Skyler and Barry 1990) of wind tunnel 
tests conducted principally at the University of California, Davis, California. Facilities to 
conduct wind tunnel tests must include: . 


1. A wind tunnel of sufficient diameter to allow horizontal expansion of the 
spray plume and of sufficient length and design to minimize turbulent air flow; 


2. Instrumentation such as a laser to measure and count particles accurately; 

3. A mechanical system for air flows up to 70 m/s; 

4. Devices to position, hold, and move the atomizers; 

5. An exhaust system to collect the exiting spray; and 

6. A pressure metering system to transport the spray to and through the atomizer. 

Other facilities known to exist include wind tunnels at Cranfield Institute of 
Technology, England; University of New Brunswick, Canada; USDA A gricultural 


Research Service, College Station, Texas; New Mexico State University, Las Cruces, 
New Mexico; and SpraySearch Daratech, Australia. 


The purpose of the present paper is to consider the physical processes that 
generate the spray drop size distribution, and to deduce simple relationships between 
nondimensional parameters that control the breakup of a jet of liquid released through a 
nozzle mounted on a spray aircraft. This work extends the summary discussion on non- 
Newtonian fluids previously presented in Teske, Bilanin and Barry (1993). 


Jet breakup atomizes the spray material into its drop size distribution. While 
many of the fundamental principles involved are well known, the results (for the most 
part) are understood only for Newtonian fluids. It is the application to non-Newtonian 
fluids that has received scant quantitative attention, and will be addressed in a limited 
sense here. 


Several researchers have investigated spray distributions generated by water- 
based and non-Newtonian fluids (Yates, Akesson and Cowden 1984; Bouse, Carlton and 
Jank 1988; Bouse, Kirk and Bode 1990; and Bouse 1991 to name a few), and have also 
tried to correlate their data (Picot, van Vliet and Payne 1989; Womac, Williford and 
Hanks 1990; and Hewitt 1993). Many other studies have been conducted (such as Esterly 
et al. 1993), but their results are unpublished or proprietary. 


Generic results from wind tunnel tests indicate the following: 


1. Shear across hydraulic nozzles caused by nozzle orientation, and shear across 
the rotating basket cages of rotary nozzles caused by rotation, are major factors that 
produce small drops; 


2. An increase in surface tension increases drop size; 
3. Viscosity is a minor influence on atomization; 


4. Rotary atomizers can become overloaded with liquid volume, resulting in 
larger drop size atomization even when the rotation rate is held constant; 


5. Most nozzles produce a large number of small drops; however, these drops 
generally represent less than one percent of the total volume in the:spray; and 


6. Slight changes in some chemical, physical and biological properties of a tank 
mix can significantly alter the atomization process. 


It is this last observation, and a concern for the smallest drop sizes (those that can 
infect the insect but also drift off-target), that underscore the problem. The need for wind 
tunnel evaluation of tank mixes and spray devices will continue to increase as new 
formulations of biological and other more environmentally acceptable pesticides are 
developed. Production and application costs dictate a need to improve application 
efficiency, efficacy, and total accountantacy and environmental fate of released material. 
Low volume applications require highly efficient atomization. Thus, there will be a high 
demand to characterize sprays in wind tunnels or through data compression. Precise 
atomization data will also be required for computer model simulations of aerial 
application. 


The difficulty in discerning the atomization of current agricultural materials is that 
the basic spray material may rarely be treated as water, and if it is, any additives to the 
tank mix quickly render the spray non-Newtonian in nature. Therefore, the real problem 
of interest here is non-Newtonian. 


2. A Non-Newtonian Primer 


_Rheology, the study of deformation and flow of matter, has traditionally grounded 
itself in Newtonian fluids. Only recently has research expanded to include non- 
Newtonian effects. The addition of polymers to pesticides, and the formulation of spray 
materials that do not use water or oil carriers, have brought home the need to investigate 
this extremely complicated subject further. In addition, the available literature suggests 
that non-Newtonian fluids cannot be approximated with Newtonian equations; small 
corrections to Newtonian behavior do not yield non-Newtonian behavior. 


Few recent publications exist that deal specifically with non-Newtonian effects. 
Lee (1979) discusses some of the more intriguing properties of non-Newtonian fluids. 
Knell (1989) hypothesizes a representation for non-Newtonian viscosity based on 
Newtonian fluids, while Mansour and Chigier (1993) attempt to measure extensional 
viscosity, Shavit and Chigier (1993) examine dynamic surface tension, and Kitamura and 
Takahashi (1988) and Dooher et al. (1988) study non-Newtonian atomization effects. 
Chhabra (1993) examines the behavior of particles in non-Newtonian fluids; his Chapter 
2 summarizes non-Newtonian fluid behavior, and contains an extensive bibliography. 


An elastic fluid may be said to possess "perfect memory”: as long as the material 
is not stressed to yield, it will return to its original state when the disturbing force is 
removed. Viscoelastic fluids (non-Newtonian fluids) possess what is known as "partial 
memory”: these fluids do not necessarily return to their original state after a disturbance 
is applied, and then removed. 


For the most part the laws of Newtonian physics are suspect for non-Newtonian 
fluids. Among the more dramatic effects these fluids exhibit are the following (Lodge 
1964 and Bird, Armstrong and Hassager 1977 with pictures): 


1. Recoil: when a non-Newtonian fluid is poured from one beaker to another, 
and the stream is cut by a pair of scissors, the upper portion of fluid already out of the 
upper beaker will retreat back into the beaker (all Newtonian fluids will fall into the 
lower beaker because of gravity). 


2. Rod-Climbing: a non-Newtonian fluid will climb a rotating rod, contrary to 
centrifugal force effects (Newtonian fluids form a depression around rotating rods). 


3. Die-Swell: a non-Newtonian fluid leaving a cylindrical tube through an orifice 
will expand to a diameter considerably larger than the tube diameter (Newtonian fluids 
are subject to the venturi effect, reaching a diameter smaller than the orifice diameter 
before expanding radially). 


4. Tubeless Siphon: when a siphon is lifted out of a beaker of non-Newtonian 
fluid, the fluid continues to be siphoned (Newtonian fluids would stop flowing). 


5. Drag Reduction: most commonly, the addition of a polymer at only 10 to 50 
parts per million may reduce the drag of a fluid by up to 50 percent over its Newtonian 
value. 


With spray materials now, as likely as not, possessing non-Newtonian behavior, it 
is clearly important to analyze the effects of these fluids on the generation of spray drop 
size distributions, and what effects they may cause beyond Newtonian behavior. 


3. Constitutive Equation 


Newtonian fluids (such as water) behave under deformation in the well-known 
linear manner: 


oU 
U— = T (1) 
or 


where i is the absolute viscosity, dU/dr is the rate of strain, and t is the resulting fluid 
stress. When the fluid is non-Newtonian, Eq. (1) no longer applies. Rather, the 
assumption is made (Goldin et al. 1969) that the rheological properties of any 
incompressible isotropic viscoelastic material in the limit of very small deformations 
follows the constitutive equation developed by Biot (1954) and Bland (1960), written for 
a fluid as: 


c 


au de NN (t-t) ]de 
U—=Ee+n—+ c exp| 2 |X ae (2) 
or >» , Th ot 


where U is the fluid velocity, r is a length, € is the linear strain, E is the modulus of 
elasticity, T| is the long-term viscous flow effect, t is ttme, N is the number of degrees 
of freedom in the fluid, c, are amplitude coefficients, and t, are time constants. The 


quantity de/ot denotes the rate of strain; E, 1, c, and t, are real and positive. c, 


and T, are present because of the non-Newtonian nature of the material. These 


quantities must be a part of any analysis or interpretation of non-Newtonian fluid 
behavior. 


4. Dimensional Analysis 


_ Dimensional analysis is used here to infer the nondimensional groups that control 
the diameters Do.1» Do.5, and Do.9 describing the structure of a drop size distribution 
(the diameters below which 10, 50, and 90 percent of the spray material is found, 
respectively). The principles of dimensional analysis are explained in any first book on 


engineering (for example, Kuethe and Schetzer 1959), and are reviewed here before 
proceeding. 


_ If an algebraic equation expresses a relationship among physical quantities, the 
equation can have meaning only if the terms involved are dimensionally identical. The 
requirement of dimensional homogeneity in physical equations is useful in determining 
possible combinations in which variables can occur. If it is required that all terms in an 
equation be pure numbers (a general assumption), then the variables involved in the 
equation may occur only in combinations that have zero resultant dimensions. Any 
physical equation can in fact be expressed in terms of dimensionless combinations of the 
variables that influence the behavior of the equation. The formal statement of this fact is 


embodied in the Buckingham II theorem; simply stated, any function of N variables: 


Fon (Py, P>,P3,P4,,...,Pn) =0 (3) 
may be expressed in terms of N-K II products: 
fen (11; , My, 13,y,..., In.) =0 (4) 


where each II]; parameter is a dimensionless combination of an arbitrarily selected set of 
K variables, and one other. K_ is equal to the number of fundamental dimensions 
required to describe the variables P;. Here all quantities P; may be expressed in terms 
of mass M, length L, and time T, and therefore K =3. The arbitrarily selected set of 
K variables may contain any of the quantities P; with the restriction that the K set 
itself not form a dimensionless combination. 


This approach cannot determine the values for the constants or coefficients in the 
functional relationships that are developed. Experimental (in this case wind tunnel) data 
must be analyzed in order to deduce the values for the constants or coefficients in these 
relationships. 


In the present problem we consider a simple jet issuing from a nozzle oriented at 
an angle into the free stream air flow (Figure 1). The variables that control jet breakup 
into drops are summarized as follows: 


VARIABLE DESCRIPTION DIMENSIONS 


D j jet or nozzle diameter L 

p fluid density ML-3 

Mm fluid viscosity ML-!T-! 
Ch stress relaxation amplitude ML-! 1-2 
Th fluid relaxation time scale T 

U; fluid jet velocity ila 

(a ee air density ML-3 
(ie air viscosity ML"!T-1 
(Whe air velocity LT! 

o surface tension MT-2 

B angle between airflow and jet - 

Q rotation rate (rotary atomizers) az! 

Dok scaling diameter L 


Jet characteristics are denoted by the j subscript. The angle B is measured such 
that B =O places the jet co-flowing with the air stream. Only fluid mechanics is 
investigated here, including the non-Newtonian effects present with c,, and T, , as given 


in Eq. (2). Any heat transfer effects because of temperature change must be reflected in 
changes in the fundamental variables, and not in a temperature variable itself 
(temperature is not typically measured when drop size distributions are determined). 


Dimensional analysis suggest that these 13 variables can form ten arbitrary 
parameters. One of them is the desired scaling parameter: 


AO = fen ( (5) 
i 


where k takes on the values of 1, 5,or 9. The rest of the functional forms may be 
obtained from previous experience and/or extended from the analysis of Lee (1979) to jet 
flow: 


IT EXPRESSION INTERPRETATION 








U; 
Jy ere 
TA velocity ratio 
B exit angle 
QD; 
II, rotation rate ratio (rotary atomizers) 
J 
= density ratio 
= viscosity ratio 
pU;D; 
hs Re —_ Reynolds number (aerodynamic force / viscous force) 
es | 
= We Weber number (aerodynamic force / surface tension force) 
D: 
a Il, non-Newtonian stress / Newtonian stress 
J 
TU; 
~ De Deborah number (liquid relaxation time / flow time) 
J 


The U; velocities present in the above expressions are interchangeable with U,, 
through the velocity ratio. The effects of hydrostatic head (through the introduction of 
gravity and the Froude number) and of compressibility (through the speed of sound and 
the Mach number) are assumed negligible in the present analysis. For water the density 
ratio p / p., = 815.3 ; specific gravities of spray material do not vary significantly from a 
value of 1.0 , implying that the density ratio need not concern us here. Likewise, the 
viscosity ratio u / U,, = 184.2 and is also not examined here. 


The Reynolds and Weber numbers are familiar to atomization processes (see 


throughout Lefebvre 1989). The unnamed parameter I1, has received no attention as far 
as our literature search could indicate, yet it is the parameter that reflects the effects of 
fluid stress. The Deborah number, on the other hand, tends to reflects the time-scale 
response of non-Newtonian fluids (Reiner 1964 and Metzner, White and Denn 1966). At 
low Deborah numbers the material exhibits a fluid-like behavior, while at high Deborah 


numbers the material exhibits a solid-like behavior. An example usually cited in the 
literature (Lee 1979) is to move a fist slowly into containers of Newtonian and non- 
Newtonian liquids -- both exhibit the same effect of allowing the hand to enter. If, 
however, the fist is thrust quickly into the Newtonian fluid, splashing results. In the non- 
Newtonian fluid the surface of the material responds much like a solid surface -- the hand 
would not enter the material and splashing would not occur. For a Deborah number 
greater than about 0.05 , spray droplet sizes can be ten percent larger than with 
Newtonian fluids. 


Any expression constructed from these parameters must have a basis in 
Newtonian fluids. Thus, it may be seen (also from Eq. (2)) that in the Newtonian limit, 


Cc, =0 and tT, =0. 


Non-Newtonian terminology sometimes includes the expressions dynamic surface 
tension and elongational viscosity (briefly discussed in Lodge 1964). Purely as an 
exercise in recovering these types of variables, it may be seen that a product of the two 
non-Newtonian parameters gives: 


Cn Ty 
rv 


which would enable the numerator of this parameter to be interpreted as an additional 
viscosity parameter (with units of ML” Iy-1 ). Also, the product: 


D: 


ay 


has the units of MT~2 , consistent with surface tension. This parameter could then be 
interpreted as an additional surface tension. What may happen, unfortunately, is that c, 


and %,, are related in complicated ways to dynamic surface tension and elongational 


viscosity. Alternately, the surface tension o and the fluid viscosity 1 may be re- 
interpreted in the Reynolds and Weber number regimes present here. Mansour and 
Chigier (1993) suggest that fluid properties should only be measured under conditions 
consistent with atomization, but readily acknowledge that suitable instruments are still 
being developed to measure the required physical properties (which in turn have not been 
decided upon). 





Figure 1. Schematic of jet breakup. 
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5. Available Data Sets 


Wind tunnel data must be available to determine the constants and coefficients 
forming the anticipated drop size equation: 


DOaefen Og Perle Codey Pe Jaap cUr ape) (6) 
Dj U7 (Uj eo Uy 


oso 


The largest available data set is a collection of all USDA Forest Service data 
(Skyler and Barry 1990), but other data sets have been published by Bouse and his 
colleagues, and by Picot and his colleagues. Other researchers (references cited earlier) 
did not vary test conditions sufficiently, and their data would be difficult to interpret here. 
The overriding requirement is that enough distinct tests be run to enable the 
nondimensional parameters to appear in a variety of values to render the definition of fen 
in Eq. (6) possible. The number of wind tunnel tests must be large enough to validate and 
test the constants and coefficients that generate the functional relationships. Suggestions 
as to the minimum number of tests needed are developed later. 


The tests for which unique data are at least possible from the references available 
(the extensive Spray Drift Task Force data set is not available at this time) are 
summarized in Table 1. The flat fan and hollow cone nozzles will be examined in detail 
in the next section, followed by Beecomist and Micronair AU5000 rotary atomizers. 


Ti 


Table 1. Candidate Wind Tunnel Data Sets 


NOZZLE TYPE 


From Skyler and Barry (1990): 


Beecomist 


Micronair AU5000 


8001 
8003 
8004 


8006 
8010 
8020 
D2-23 
D2-25 
D4-45 
D4-46 
D8-45 
D8-46 
D8 Jet 
RD-7 
RD-10 


MATERIAL 


Dipel 8L neat 
Dipel 8L 1:1 
Foray 48B neat 
T Biocontrol 
Water 

Dipel 6L 1:1 
Foray 48B neat 
Thuricide 32LV 1:1 
Water 

Water 

Foray 48B neat 
Foray 48B neat 
Water 

Water 

Water 

Water 

Water 

Water 

Water 

Water 

Water 

Water 

Water 

Water 

Water 


From Bouse, Carlton and Jank (1988): 


D6-45 


D6-46 


Water with Nalco-Trol 
Water with Nalco-Trol I 
Water with Sta-Put 
Water with Drifgon 
Water with 38-F 

Water with Wind-Fall 
Water with Nalco-Trol 
Water with Nalco-Trol I 
Water with Sta-Put 
Water with Drifgon 
Water with 38-F 

Water with Wind-Fall 


From Picot, van Vliet and Payne (1989): 


Micronair AU5000 


Water with Isopropanol 
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NUMBER OF TESTS 
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6. Water: Flat Fan and Hollow Cone Nozzles 


The most obvious first step is to investigate the flow of a Newtonian fluid, water, 
through the flat fan and hollow cone nozzles. The available data sets (from Skyler and 
Barry 1990) yield the wind tunnel results summarized in Table 2. Note that the data for 
the 8006 (Table 1) are not included; we will measure the performance of the technique by 
comparing its predictions with the 8006 data. Nozzle diameters are obtained from the 
latest Spraying Systems Company Industrial Spray Products Catalog (No. 51). 


An initial examination of the data revealed that much better correlation with data 
was possible when the velocity present in the Reynolds and Weber numbers is the 
relative velocity: 


U., - Uj 


without a cosine correction for angle B (the cosine correction makes things much 
worse). Therefore, for the flat fan and hollow cone nozzles, we define: 


pCU,-Uj)D; 
Re = (7) 
vy 
p(U,.-Uj)*Dj 
We =—_____ (8) 
o 
The assumed functional relationship is: 
aes = A, + By Rea Wed (9) 
J 


where k=1, 5,o0r 9, and the dependence on Re and We is universal across all drop 
size distributions. This further assumption (the constancy of a and b ) permits us to use 
all three sets of diameters recovered from the wind tunnel tests to develop a least squares 
algorithm and minimize the standard deviations between the data and the predictions of 
Eq. (9), while also maximizing the correlation coefficient (Guttman, Wilks and Hunter 
1971): 


2_ 2 (P-P)(D-D) 
Sp OD 


R (10) 
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where D and P are the data and predicted values, respectively, Op and Op are the 


standard deviations of the data and predicted values, respectively, and an overbar denotes 
an average value. 


Exploration of all reasonable values of the powers on the Reynolds and Weber 
numbers suggest that a broad region exists where the correlation coefficient, as computed 
from Eq. (10), remains fairly constant. A quantitative test could not be constructed to 
discern the "best" values for a and b. Therefore, it was decided to fit the available data 
to an a= -1.0 power law on Reynolds number. This assumption suggests that 
agricultural spray atomization behavior is consistent with other atomization applications 
(Lefebvre 1989) which recover a parametric variation with the Ohnesorge number: 


_ We9-5 
Oh = We" (11) 


In addition, it was decided that the power law on Weber number should be the 
same value irregardless of nozzle orientation. This process leads to the representations: 


FLAT FANS AT 8 = 0 DEGREES (see Figure 1) 


70.1 = 0.0198 + 61.84 Ree (12) 
‘| € 

70.5 = 0.1101 + 80.44 [weed £2) (13) 
‘| € 

P09 - 0.2290 + 99.61 we. A2 (14) 
Dj e 


with R2 = 0.849, 0.645 and 0.529 respectively. These results are not as well-correlated 
as with the other nozzle angles. 


FLAT FANS AT B = 90 DEGREES 


DO.1 — -9,0217 + 69.70 ELE (15) 
Dj Re 
v0.5 = -0.0318 + 129.50 [wel.4? (16) 


J 
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70.9 = 0.0186 + 179.98 we. 42 (17) 
j e 


with R2 =0.955, 0.965 and 0.906 respectively. 


FLAT FANS AT £8 = 135 DEGREES 


DO.1 — 9.0201 + 65.31 [wed-*2 (18) 


DO.5 — .0,0491 + 129.08 cues (19) 
Dj Re 





DO.9 ~ -0,0726 + 189.50 cag (20) 
D; Re 


with R2 = 0.934, 0.978 and 0.974 respectively. 


HOLLOW CONES AT 8 = 0 DEGREES 


re = 0.0143 + 82.12 wet 28 (21) 
DO.5 = 0,0352 + 139.00 [we 38 (22) 
Dj Re 
70.9 = 0.0731 + 207.20 jweed.38 (23) 
J 


with R2 = 0.874, 0.757 and 0.598 respectively. These results are not as well-correlated 
as with the other nozzle angle. 


S 


HOLLOW CONES AT B = 90 DEGREES 


PO.1 = 9.0015 + 83.58 we 38 (24) 
j 
PO.5 = 0.0011 + 149.37 | We0-38 (25) 
D; Re 
20.9 = 0.0025 + 228.71 wed.38) (26) 
j 


with R2 = 0.891, 0.918 and 0.857 respectively. 


These correlations are plotted in Figures 2 - 6. 


The above results are not consistent with other researchers who have attempted 
data compression (Picot, van Vliet and Payne 1989 and Hewitt 1993). Their difficulty 
was in trying to correlate the effects of non-Newtonian fluids immediately, without first 


developing the Newtonian effect. 


As an example of the accuracy of the above equations, we next predict drop size 


distribution with the Spraying Systems 8006 nozzle at 6B = 135 deg (Table 1), which was 
not used to generate the constants and coefficients appearing in Eqs. (18) - (20). This 
example gives: 


VARIABLE PREDICTION DATA DIFFERENCE ERROR 


Do 1 (um) 89.29 100.34 11.05 11.0 % 
Do 5 (um) 161.69 175.40 13471 7.8 % 
Do 9 (um) 236.52 257.73 21.21 8.2 % 


Overall, the accuracy of the curve-fit equations is estimated to be 14 percent. 
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NOZZLE 


8004 


8006 


8010 


8020 


ANGLE 


AIR SPEED 


bil 


49.2 
53.6 


22.4 


67.1 
22.4 


60.4 
67.1 
22.4 


49.2 
67.1 


* 22.4 


67.1 
22.4 
44.7 
67.1 
22.4 
44.7 
67.1 
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Table 2a. Wind Tunnel Data for Water: Flat Fan Nozzles 


Do.5 
(4m) 


204.77 
212.10 
197.51 
2ozeey 
181.56 
139.01 
218.67 
180.54 
141.97 


362.55 
333.64 
258.51 
SEERA 
242.17 
169.72 
296.65 
201.70 
ey hehe’ 


319.84 
175.40 


542.34 
518.31 
290.82 
413.98 
2595) 
221.11 
182.37 
343.93 
Z10o 
212.05 
173.49 


686.11 
609.67 
332.18 
478.49 
294.51 
211.80 
394.48 
236.83 
184.56 


1030.60 
766.56 
571.14 
697.76 
408.01 
B72 15 
Zales 
537.49 
316.41 
318.42 
237.94 


1001.48 
996.89 
620.58 
797.44 
483.75 
425.30 
616.75 
358.26 
262.42 


NOZZLE 
iveys 


D2-23 


D2-25 


D4-45 


D4-46 


D8-45 


D8-46 


Table 2b. Wind Tunnel Data for Water: Hollow Cone Nozzles 


ANGLE AIR SPEED 


(deg) 


(m/s) 


22.4 
44.7 
67.1 
22.4 
44.7 
67.1 


22.4 
44.7 
67.1 
22.4 
44.7 
O12 


22.4 
44.7 
67.1 
22.4 
44.7 
67.1 


22.4 
44.7 
67.1 
22.4 
44.7 
67.1 


22.4 
44.7 
67.1 
22.4 
44.7 
67.1 


22.4 
44.7 
67.1 
22.4 
44.7 
67.1 
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7. Water: Rotary Atomizers 


The next step is to investigate the flow of a Newtonian fluid, water, through rotary 
atomizers. Again, the available data sets (from Skyler and Barry 1990) yield the wind 
tunnel results summarized in Table 3. Note that the data contains both Micronair and 
Beecomist atomizers. Details of the Micronair AU5000 are obtained from Picot, van 
Vliet and Payne (1989), and for the Beecomist, from P. G. Phillips (private 
communication). 


Both Picot, van Vliet and Payne (1989) and Parkin and Siddiqui (1990) correlated 
rotary atomizer data and found their best results when the characteristic size of the 
atomizer is taken to be D, , the diameter of the cylindrical mesh shell. Also, Parkin and 
Siddiqui (1990) found that the rotation rate best correlates with their data as the inverse 
square root. This information will be used here to extend the results presented above to 
rotary atomizer behavior. 


Therefore, for rotary atomizers, we define: 


p(U,,-Uj)De 
Re = ———_ 21) 
tv 
p(U,,-Uj)*D¢ 
We = ——_—_—— (28) 
oO 
The assumed functional relationship is: 
as = Ay + By Re? Wed / 11,0 (29) 
Cc 


where k=1, 5, or 9, and the dependence on Re and We is again universal across all 
drop size distributions and rotary atomizer types. This assumption (the constancy of a 
and. b ) permits us to use all three sets of diameters recovered from the wind tunnel tests 
to develop a least squares algorithm and minimize the standard deviations between the 
data and the predictions. This process leads to the representations: 


ROTARY ATOMIZERS 
DOQ.1 — 0.00083 + 303.57 vo (30) 
De Re II," 
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DO.5 = .0,00122 + 509.60 eas (31) 


Do Re 15 
20.9 = -0.00167 + 721.03 Py (32) 
Cc Re I,’ 


with R2 = 0.813, 0.897 and 0.861 respectively. The correlation is plotted in Figure 7. 
Parkin and Siddiqui (1990) apply a least-squares algorithm to their data to arrive 
at the power-law dependences: 


U..2° 18 


p, 0.70 


while the above results give: 


OnUpe 
D,7!.27 


Our analysis suggests an attempt by the nondimensional parameters to "remove" air speed 
as a variable in rotary atomizers. 


(A) 


Table 3. Wind Tunnel Data for Water: Rotary Atomizers 


NOZZLE ROTATION AIR SPEED Do1 Do5 Doo 
TYPE (tpm) (m/s) (um) (Um) (Um) 
AU5000 3750 22.4 106.81 248.73 284.02 
9100 44.7 56.36 110.17 162.12 
11700 58.1 49.01 81.37 138.52 
4100 22.4 142.01 306.49 404.20 
8000 44.7 56.35 118.19 179.08 
10850 58.1 44.08 79.13 131.12 
4200 44.7 96.31 200.18 278.70 
7000 44.7 104.84 162.16 219.84 
11000 60.3 52.64 110.54 174.23 
Beecomist 11700 22.4 55.97 124.82 173.30 
11700 44.7 59.65 123.46 170.42 
8300 22.4 111.13 191.05 285.92 
8350 44.7 112.65 178.25 270.46 
9000 24.6 128.27 187.85 282.84 
7000 24.6 147.52 248.95 382.27 
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8. Non-Newtonian Fluids 


_ The addition of viscoelastic effects requires us to expand the functional 
relationship of Eq. (9) to read: 


PO.K =[ Ay + By Ret Web][1+C, 11,°][ 1 + Dy Ded] (33) 


Dj 


where Ay, By, a and b are defined previously (these constants and coefficients 
should not change when fluids other than water are used; the rotation rate behavior is also 
added for rotary atomizers). Because of the nature of non-Newtonian fluids, we expect 
Cy, and Dy to be positive numbers. We identify the two non-Newtonian parameters by: 


CnDj 
(Rb a 
$n _-Up ey 
_ T(U eit} ) 

e€= 2 aD ws (35) 
An alternate expression for the non-Newtonian terms could be: 
[1+E, 1,° Ded | (36) 


combining the individual expansions. Sufficient wind tunnel data will be needed to detail 
which terms to use here. 


The existing data bases do not include enough data to define quantitatively the 
new constants and coefficients in Eq. (33) above. Two difficulties arise when attempting 
to make use of the existing data sets: 


1. The viscoelastic fluid parameters c, and t, are not known, or have not been 


determined for any spray material in the data base. However, these constants could be 
worked into modifications of the constants C, and Dy , with the use of the data sets to 


"back out” their values. Unfortunately, this then leads to: 


2. There are too few unique wind tunnel tests to close the problem. For the least 
squares technique to work with some degree of confidence, there should be at least twice 
as many data points as unknown constants or coefficients. In the case of viscoelastic 


fluids through flat fan and hollow cone nozzles, even with c, and Tt, known, four 
constants must be evaluated, requiring eight wind tunnel tests. With the viscoelastic 
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physical properties not known from separate benchtop tests, an additional constant must 
be determined, requiring two additional wind tunnel tests. 


Physical properties of the spray material must be determined, then additional wind 
tunnel tests must be performed, before we can quantify the behavior of non-Newtonian 
fluids in this application. 


Picot, van Vliet and Payne (1989) contains seven tests with the Micronair 
AU5000 for a non-Newtonian fluid, water with nine percent isopropanol (Table 4). This 
paper also includes measurements of the surface tension and viscosity. If we use the 
alternate expression (Eq. (36) above) and expand the non-Newtonian terms, we obtain the 
expression: 





WPL lee (as S y (37) 
2 a ircomescp] De 
with: 

F, = Ex Cac td (38) 


This approach contains three constants: the powers c and d and the coefficient F,. 


However, all of the data is collected at the same air speed, which suggests that the terms 
in parentheses in Eq. (37) will change very little in a least-squares analysis. This effect 
further suggests that the best possible examination of this data would be to assume that 
the non-Newtonian terms may be approximated by the behavior: 


[1+G,] | (39) 


essentially a constant multiplying the Newtonian result. This approach leads to the 
values of: 


Gy = 5.85 for Do.1 
Gy = 3.90 for Do5 (40) 
Gy = 3.39 for Do.9 

with R? = 0.754, 0.834, and 0.788 respectively. These results are encouraging, but 


indicate the need for further study and laboratory or wind tunnel data. Equation (40) 
would suggest that the non-Newtonian effect is fairly significant. 
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Table 4. Wind Tunnel Data for Water/Isopropanol: Rotary Atomizers 


NOZZLE ROTATION AIRSPEED Do; Do.5 Do.9 
TYPE (tpm) (m/s) (Lum) (um) (14m) 
AU5000 7920 44.7 42.2 106.0 182.0 
11000 44.7 41.1 94.0 175.0 

3660 44.7 57.7 165.0 261.0 

7920 44.7 43.9 126.0 216.0 

11000 44.7 47.8 115.0 195.0 

3660 44.7 71.7 161.0 284.0 


7920 44.7 49.2 136.0 239.0 
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9. Drop Size Distributions 


With the Dp. 1, Do.5 and Dog drop diameters in hand, we may then turn to 


the construction of the entire drop size distribution profile, to obtain an indication of the 
accuracy we have implied in the dimensional analysis approach. As an example here, we 


will use the 8006 flat fan nozzle result (at an exit angle of B = 135 deg ), as given in the 
Skyler and Barry (1990) data base (and plotted in Figure 8), whose cumulative volume 
fraction CVF is plotted in Figure 9. For this example, Table 2 indicates that: 


Do,1 = 100.34 um 
Do.5 = 175.40 um 
Dog = 257.73 bm 


Several techniques have been developed to curve fit the wind tunnel drop size 
distribution. The most widely used expression for displaying drop size distribution is the 
logarithmic normal representation (Rosin and Rammier 1933): 


1 - CVF = exp[ -(D/X)4]| (41) 


where X and q are constants, and D is the drop diameter. Traditionally, X is taken 
as the drop diameter Do 632 . and q is the slope of the least squares curve fit through 


the data in transform space. For the 8006 data, transform space may be plotted as shown 
in Figure 10 (with the parameters Do 637 = 198.47 um and q = 3.1316 ). The resulting 


drop size distribution comparison is given in Figure 11 ( R2 = 0.978 ), and is seen to miss 
at the maximum volume fractions (near Do 5 ) and be high at the larger drop diameters. 


A more complicated technique is the upper-limit function (Mugele and Evans 
1957), requiring the determination of three constants correlating CVF to an error function 
behavior. Goering and Smith (1978) give approximate formula for generating these three 
parameters: 


Do.5(D9.1+D99)-2D91D 
ee Oe eee ta 0,9 . (42) 
Pm Do.5 - Do.1 Do.9 / Do,5 
A =Dm-50.5 (43) 
Do.5 


Sat Reet 0204 2 
ay oy SAND of 
As Dm -Do.9 (44) 
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to give 


CVF = 0.5 (1+erf(5z) ) (45) 
=log (Bi 
2= log (AOL) (46) 


and erf is the commonly-defined error function. For the 8006 data, D,, = 396.40 um, 


A = 1.260 , and 5 = 1.066. The correlation to CVF is shown in Figure 12, and the drop 


size distribution comparison is given in Figure 13 ( R2 = 0.987 ). The curvefit is closer to 
the data than the Rosin-Rammler representation of Figure 11. This approach has recently 
been extended to four constants for better accuracy (Patel and Gaidos 1991). 


The technique preferred here (Teske, Skyler and Barry 1991; Teske and Barry 
1992; and Teske 1992) was developed by Simmons (1977) from an idea he attributed to 
Tate and Marshall (1953). In this technique CVF is plotted on a normal probability 
distribution scale as a function of the square root of the drop diameter. The results were 
generally so startling that the authors felt compelled to dismiss divine intervention and 
intimate that the correlation was purely empirical and should not be used as a basis for 
theoretical interpretation. 


If the drop diameters are normalized by Do 5, and a least-squares fit is invoked 
in root/normal space, with an iteration technique employed to maximize the area under 


the volume fraction (Esterly 1991), Figures 14 and 15 ( R2 = 0.989 ) result. In this case 
only two parameters are needed to describe the distribution: Do 5 and the slope a of the 


straight line in Figure 14. Here Do 5 = 176.47 um and a=0,1738 , and only the tail of 
the distribution is predicted incorrectly. 


If we take the 8006 data and realize that: 


A / Doo . 4 | Do. 
ae D0 Sas” (47) 
Pr(0.9) -Pr(0.1) 
where Pr() is the normal probability distribution function (forwards and backwards in 
equations 26.2.18 and 26.2.23, respectively, of Abramowitz and Stegun 1968), for which 
Pr(0.9) = -Pr(0.1) = 1.282. In this application the slope parameter becomes a= 0.1778, 
an especially close fit to the slope in Figure 14. 


Alternately, the relative span RS , defined by: 


oe. 


— 20.9 - Doi (48) 
Do.5 ; 


may be supplied. Through manipulation of the root/normal function, Eq. (47) may be 
replaced by the solution to the equation: 


eo Roe 
a= TPr05) 7 0-195 RS (49) 


For this case, RS = 0.897 from Eq. (48), and a=0.1749 from Ret (49), an even better 
fit to the slope in Figure 14. 


oe ate use the root/normal approach on the predicted diameter values (through Eqs. 
18) - (20)): 


Do.1= 89.29 um 
Do.5 = 161.69 um 
Dog = 236.52 pm 


which give a=0.1776 from Eq. (49), the predicted drop size distribution that results is 


shown in Figure 16 ( R2 = 0. 932). Here it may be seen that the prediction closely 
matches the actual wind tunnel data. 
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10. Procedures 


It would seem prudent now to review the procedures that should be followed to 
correlate drop size distributions across available wind tunnel data sets: 


1. For the nozzle types for which a correlation already exists, the researcher 
should first compare the Newtonian predictions with the Newtonian equations (for 
example, Eqs. (18)-(20)) for the spray material in question. If the predicted drop sizes are 
within ten percent of the data, the assumption of Newtonian behavior (as represented in 
the drop size distribution) may be a good one. 


2. If not, then the viscoelastic physical properties must be measured, and the 
non-Newtonian equations used (as developed, for example, in Eq. (33)). 


3. If the nozzle has not been parameterized, then water should first be wind- 
tunnel tested, to develop the Newtonian functional relationship and evaluate the constants 
_and coefficients. Then the spray material in question should be wind-tunnel tested, and a 
comparison of data made with the Newtonian equations just developed. If the predicted 
drop sizes are within ten percent of the data, the spray material may be taken as 
Newtonian. Otherwise, additional wind tunnel tests will be needed to develop the non- 
Newtonian functional relationship and evaluate the additional constants and coefficients. 
Physical properties will be needed to close the problem accurately. 


4. It is important that accurate physical property measurements and drop size 
distributions be obtained. Otherwise, the correlation coefficient for data compression, 
reflecting all of the errors built into the data, will be low. 


5. It is important that sufficient wind tunnel tests be obtained to make the least 
square approach discussed here work well. For all the data sets in the USDA Forest 
Service data base, and all other data found in the literature search, few tests (beyond 
water) were performed in the systematic way necessary to invoke the techniques 
discussed here, and enable an accurate prediction of the drop size distribution beyond the 
wind tunnel test conditions. Limiting wind tunnel tests must be included in any study to 
cover the entire anticipated range of the application for the spray material. Clearly, a 
minimum of eight distinct wind tunnel tests are needed to extract the four Newtonian 
constants (running water) and eight additional tests are needed to extract the four non- 
Newtonian constants. The cost of this testing must be weighed against the use of the 
spray material in the field application: it makes no sense to develop the non-Newtonian 
equation if the spray material is always applied in a specified concentration at a specified 
nozzle angle by an aircraft spraying at a specified speed. 
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11. Conclusions and Recommendations 


_ This study has reviewed the state of Newtonian and non-Newtonian spray 
material characteristics, and attempted data compression with available wind tunnel data. 
Starting with first principles of dimensional analysis, we have developed consistent 
nondimensional parameters that group into well-accepted sets (Reynolds and Weber 
numbers, for example). We have correlated the water runs from the USDA Forest 
Service drop size distribution data base with a simple expression that recovers the details 
of the spray for conditions that were not included in the equation development. We have 
ee upon the form of the non-Newtonian terms for the atomization of viscoelastic 

uids. 


What is now needed is a two-fold process: 


1. A systematic laboratory examination of non-Newtonian fluids (and their 
Newtonian counterparts) should be undertaken. Benchtop testing should be considered to 
better understand the behavior of viscoelastic fluids, using small quantities of 
commercially available non-Newtonian spray materials. The development of their 
physical properties, particularly the proposed stress relaxation amplitude and fluid 
relaxation time scale constant, must be better understood. 


2. A consistent set of wind tunnel tests must be performed for anticipated spray 
materials and field application conditions. These tests must retain the high Reynolds and 
Weber number regimes present in the field spray applications, and recover the smaller 
drop size categories (down to 10 micrometers) because these are the drops most likely to 
drift, and they must therefore be quantified (through Dg 1 ) as well as possible. The 


guidelines in this paper should be followed to quantify the number of tests to conduct to 
recover sufficient data to quantify the behavior of the tested spray material. 


It is important in both recommendations that accuracy be a high priority. Any 
experimental deviation will carry through to the final equations developed. Also, it 
would be hoped that the extensive data base developed by the Spray Drift Task Force 
would be available for extending the analysis in this paper, and quantifying the effects 
suggested here. 
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